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ABSTRACT 

We present i?-niatrix calculations of photoabsorption and photoionization 
cross sections across the K edge of the Li-like to Ca-like ions stages of Ni. Level- 
resolved, Breit-Pauli calculations were performed for the Li-like to Na-like stages. 
Term-resolved calculations, which include the mass- velocity and Darwin relativis- 
tic corrections, were performed for the Mg-like to Ca-like ion stages. This data 
set is extended up to Fe-like Ni using the distorted wave approximation as im- 
plemented by AUTOSTRUCTURE. The i?-matrix calculations include the effects 
of radiative and Auger dampings by means of an optical potential. The damping 
processes affect the absorption resonances converging to the K thresholds causing 
them to display symmetric profiles of constant width that smear the otherwise 
sharp edge at the K-shell photoionization threshold. These data are important 
for the modeling of features found in photoionized plasmas. 

Subject headings: atomic data - atomic processes - line formation - X-rays: 
spectroscopy 
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Introduction 



The K lines and edges of iron and nickel are important diagnostics in X-ray astronomy. 
They are in a relatively unconfused part of the X-ray spectrum (~ 6-10 keV), and they 
are formed over a wide range of physical conditions. This can range from cold (by X-ray 
standards) neutral gas at temperatures < 10*^ K, to highly ionized gas near 10^ K. Their 
energy and shape give information about the ionization and kinematics in the gas. The 
strengths of emission features, such as K fluorescence lines, constrain the amount of emitting 
gas and the elemental abundances, while the strengths of absorption edges constrain the 
column densities along the line of sight to the background X-ray source. Iron K lines and 
edges have been studied by many existing and past X-ray astronomy instruments. They have 
provided insights into d iverse topics such as ele mental abundances in clusters of galaxies 



(ISer 



emitsos et al. 



19771 1 , the spin of blac k holes (iBrenneman et al. 



20111 ). and the geometry of stellar flares ( jTesta et al. 



2011 : 



McChntock et al. 



20081 ). Nickel features have not yet 



been used to the same extent, owing to the fact that they are weaker by a factor ~10 from 
the typical Ni/Fe abundance ratio, and to the fact that many detectors have collecting 
areas which decrease rapidly above 7 ke V. Nonetheless, Ni c kel has been detected in systems 



ranging fror n black hole X-ray binaries ( iMiller et al. 



2006 



of galaxies ( 



remnants ( JTamagawa 



Tamura et al. 



20091 ). the center of our galaxy (JKoyama 



Kallman et a. 



20091), clusters 



201 ll ). and supernova 



2010l ). Further impetus for calculations of the spectra of nickel comes 



from the Astro-H satellite, which will result in the most sensitive observations to date of 
objects at energies containing K features of nickel. This is because of both the relatively 
large sensitivity at the energy of the nickel K lines and edges (7.2 - 10.7 keV), and because 
of its spectral resolution, which will exceed that of previous experiments at energies above 
3keV. 



Over the last several years, our group has worked on providing accurate K-shell 
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atomic data for all ions of astrophysical interest. We have calculated energy levels, 
wavelengths, Einstein A-coefficients, Auger rates, and photoabsorption/io nization cross 



sections. Complete data se ts have been computec 



2003 



Palmeri et al. 



2004h. O f 



2003a 



Garcia et a. 



(^Palmeri et al. 



C sequence by 



2008a: 



b; 



Mendoza et al. 



2005h. N f 



2004: 



Garcia et al. 



Witthoeft et al. 



Hasoglu et al. 



2009 



'or all ion stages of 



Bautista et al. 



2004 



'ei 



Bautista et al 



Kallman et al. 



20091 ). as well as Ne, Mg, Si, S, Ar, and Ca 



20111 ). There has also been recent work on the 



(I2OIOI ) using similar techniques. 



The data we present in this work are important both for their direct effect on observed 
astrophysical X-ray spectra via absorption features, and also through their effect on 
ionization balance calculations, which in turn affects many of the observables from these 
elements. In particular, it is the detailed resonance and edge structure, missing in previous 
calculations, which is most important. 

There have been a large number of calculations for these systems going back to the 
1920's, but they are confined to the non-resonant background and few span the K edge 
region. For the ions considered in this work, there is good agreement between our results 
and t he background cross section both above and below the K edge characterized by the 



fits of 



Verner fc YakovlevI (119951 ). Laboratory measurements have been mainly confined to 



neutral systems for photon energies beyond the K edge . A co mprehensive assessment and 



fitting of measurements has been performed by 



Veigelel fll973h . 



There has been rece nt progress raeasur ing the res onance structure ne ar the K edge 



of Li-like sys tems: C^"*" 



Be-like C2+ fIScullv et al. 



Miiller et al 



20091 ) and B2+ flMiiller et al. 



2OIOD, as well as 



20051 ). i?-matrix calculations performed as part of those works 
are in good agreement with the measurements apart from small discrepancies in some 
resonan ce positions and he ights. With the advances of free electron lasers to produce 



X-rays (JKanter et al. 



20061 ). we expect it will be possible over the next several years to 
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perform similar experiments for the K-shell absorption by heavier species. However, such 
measurements for Ni are not yet available. 

2. Numerical methods 

The -R-ma trix method ernployed here is the same as previous K-shell calculations in 



this series, c.f. 



Bautista et al. 



appro xi mation of 



fll97lh 



Burke &: Seaton ( 



Berrington et al. 



(I1974J . 



f|20p3|). T he i?-matrix method is based on the cl ose-coupling 



1978 



19711) . w hich is solved numerically following 



Burke et al. 



19871 ). 



As in our previous calculations, the effect of r adiative and spectator Auger dec ay on 



Gorczyca fc Badnelll (11996 . 



20001). Here, 



resonances is included using the optical potential of [ 

the resonance energy with respect to the threshold acquires an imaginary component which 

is the sum of the radiative and Auger widths of the core. The radiativ e widths are compu ted 



within the i?-matrix calculations following the prescription given by 



Robicheaux et al. 



19951). The Auger w idths are taken from the relativistic Hartree-Fock calculations of 



Palmer i et al. 



(l2008bl ). These effects are referred to as radiative and Auger damping since 
they cause resonances in the photoionization cross section to become smaller in area. The 
damping also makes it necessary to calculate both the photoabsorption and photoionization 
cross sections since not every absorbed photon will result in a photoelectron. 



We use the i?-matrix computer package of 



Berrington et al. 



calculation and the asymptotic region code stgbfOdamfI^ (JGorczyca fc Badnelll 119961 ) is 



1995f) for the inner re gion 



used to compute the photoabsorption and photoionization cross sections including the 
effects of damping. 



^http://amdpp. phys.strath.ac.uk/tamoc/code. html 
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For the Li-like to Ne-like ions, our target expansion includes all states in the 
configurations: Is^ 21"^ and Is 2/^"+^ where w = for Li-like and w = 7 for Ne-like. Note 
that the target refers to the final state which has one fewer electron than the initial state; 
thus the Ne-like case having a 9-electron target. For the remaining ions, we must include 
n = 3 orbitals. The included configurations are Is^ 2/* 3/"', Is^ 2/^ 3^"^ 3d, Is^ 2/^ 3/"'+\ 
and Is 2/® 3/^"+^ where / is restricted to s and p electrons. We only include 3d excited 
states from the valence configuration to improve the atomic structure. We have tested the 
addition of 3d configurations in the L- and K-hole states for Mg-like and Ar-like Ni and 
have found the following differences. First, the KP resonances are shifted by a couple eV 
which is within the uncertainty of the calculation. This difference in the resonance positions 
decreases approaching the K edge. Secondly, the larger calculation has a few additional 
resonances, but these are small and disappear when the cross section is convolved with a 
resolution reflective of current X-ray detectors. Therefore, we are confldent in omitting the 
additional 3d conflgurations for all of the present calculations. 

We include enough continuum basis orbitals to span from threshold to at least four 
times the energy of the K edge. This choice gives reliable cross sections up to at least twice 
the energy of the K edge. Partial photoionization and total photoionization/absorption 
cross sections are calculated from all states in the ground conflguration of the parent ion. 

The close-coupling expansion is small enough for the Li-like to Na-like ions to allow 
for fully level-resolved, Breit-Pauli (bprm) calculations to be carried out. The number of 
flne-structure levels for the other ions quickly becomes too large for the computational 
resources available, so we revert to LS'-coupling for the remaining ions. The mass-velocity 
and Darwin terms are still included in these calculations which g reatly improve the t arget 



energies. Wave functions are obtained using autostructure (JBadnell 



1986 



19971 ) 



where the level/term energies have been adjusted to match those in the NIST database 
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(iRalchenko et al 



by 



Palmeri et al 



20081) w hen possible, otherwise we match our energies to those calculations 



(l2008bf ). Some levels in our target expansion are not listed by either of 
the above sources. For such valence and L-vacancy states, we use our autostructure 
energies unshifted. For the unlisted K-hole states, we shift our energies to match the 



average difference between the K-hole states appearing in both our target and 



Palmeri et al. 



(J2008bl ). This last shift is approximately 25 eV (0.3%) for the Mg-like to Ca-like ions. The 



resulting cross sections are shifted in energy so that the ionization threshold matches NIS 



his fin al shift brings the K edge thresholds within 0.5% of those given by 



Palmeri et al 



(|2008bh . 



Cross sections for Sc-like to Fe-like Ni are calculated using the distorted wave method 
as implemented by autostructure. These calculations use the same configuration 
expansion as the i?-matrix calculations with the exception of Fe-like Ni which includes the 
3d^ 4i [i = 0, 1, 2) configurations. Resonances are included using the isolated resonance 
approximation up to ra = 20. These calculations are performe d in LS coupli i ig and the 



Palmeri et al. 



toom . The 



energies are shifted so that the location of the K edge matches 
background cross section is calculated up to 13.6 keV while resonances are calculated only 
near the K edge, autostructure can construct wave functions from either the iV-electron 
target orbitals or the (A^ + l)-electron parent orbitals The difference in the cross sections 
using both sets of wave functions gives an idea of the uncertainty of this method. The 
wave functions are composed of relaxed (non-orthogonal) orbitals. Since we do not attempt 
to take radiation or Auger dampings into account, we only report total photoabsorption 
cross sections which are convolved with a Gaussian profile. The width of the Gaussian is 
AE/E = 10^^ which is represetative of currently deployed X-ray detectors. Calculations 
for Ni I and Ni II are more involved, so we leave these for future work. 



3. Results 

The cross sections reported here contain many of the same features shown in the 
results from our previous calculations. Inclusion of radiation and Auger damping broadens 
resonances near the K edge and makes them more symmetric. The damping is also 
important in distinguishing the absorption and ionization cross sections. In Figure [H 
we show these cross sections for the ^Sq ground state of Ne-like Ni. As there are no Ka 
resonances for Ne-like ions, the first resonance is K/3. The effect of damping is observed in 
both the photoionization and photoabsorption resonances. In the former, the heights of the 
resonances are greatly diminished, especially as the photon energy approaches the K edge 
due to the dominance of the spectator Auger process. For photoabsorption, both radiation 
and Auger damping broaden the resonances and make them more symmetric. 

For this simple system, one would expect a single resonance for each Rydberg n-shell, 
however, the K/3 resonance has two peaks in our results. This represents the breakdown 
of LS coupling at large Z. For \ossr-Z ions, the 2p^ ^Sq ground state goes to the Is 2s^ 
2p^ 3p ^Pi state. At high Z, the total spin is no longer a conserved quantity, making 
the ^Pi state accessible leading to the second K/3 peak. For the n > 3 resonances, LS 
coupling remains a good description for a Rydberg electron attached to a core, so only a 
single resonance is found. The appearance of the second K/3 resonance as Z increases is 



demonstrated in Fig [2 



to 28 using data from 



where we show t he K/ 3 resonance for six Ne-like ions from Z = 12 



Witthoeft et al. 



(120091 ). The separation of the K/3 resonances for 



Ne-like Ni is approximately 5 eV which should be detectable in the next generation of X-ray 
observatories, such as Astro-H. 



'1 



In Figure [3l we show the background cross sections at fixed photon energies above and 



below the K edge (at 10.9 keV and 6.8 keV, respectivel y) for all Ni ions in c 



work. For comparison, we also include the results from iVerner fc Yakovlevi (1l995l ). Both 



uded in this 
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results agree within a few percent along the whole iso- nuclear sequence. The background 
below the K edge steadily increases until the n = 2 shell is filled (Ne-like), where the 
background becomes constant. This is because the background is dominated by the L-shell 
cross section and the 2s and 2p electron wave functions are not affected much by the 
addition of n = 3 electrons. For the same reason, but for Is electrons, the background cross 
section above the K edge is relatively constant beyond He-like Ni. 

Resonances in the cross section are often quite narrow, even when Auger damping 
is included, requiring tens of thousands of energy points to resolve all the features. In 
modeling codes, this amount of data is cumbersome, so we convolve our cross sections with 
a Gaussian profile having a width of AE/ E = 10"^. An example of the convolution is 
shown in Figure H] for Be-like Ni for both photoabsorption (top panel) and photoionization 
(bottom panel). The convolved results demonstrate the strength of damping as the resonant 
contribution to the cross section for photoionization is nearly gone at the K edge. For 
photoabsorption, on the other hand, the resonances converge onto the K edge, causing the 
perceived edge in the convolved results to occur at a smaller energy; this is known as the 
smearing of the K edge in absorption spectra. 



Naharl (120051 ) performed 



As part of an effort to calculate recombination rate coefficients, 
a Breit-Pauli i?-matrix calculation for Li-like Ni with energies extending beyond the K 
edge. We are in good agreement with those results for the ground state with respect to the 
resonance positions and the location of the K edge. However, there is some disagreement 



with the background cross sectio n in the vicini ty of the edge. In Figure 



are compared with the results of 



Nahail (120051 1 and 



Slthe present results 



Verner &: Yakovlev 



(Il995h . Below the 



edge, we find fair agreement, the background of Nahar being about 10% higher than the 
present results, while above the edge, the Nahar background is nearly 30% lower. The 
disagreement above the edge disappears entirely above a photon energy of 11.5 keV. The 
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background from the present results are in excellent agreement with IVerner fc Yakovlev 



( I1995I ) at all energies, alth ough their K edge position is about 50 eV lower than that 



calculated by ourselves and 



NahaiJ (|2005|). 



We find larger differences with 



Nahan (120051 ) in the amount of resonant enhancement. 



This appears to be due to the higher energy resolution used in our calculation. Unlike the 
rest of the ions covered in this work, the Li-like system does not have a spectator Auger 
process which would significantly damp the narrow photoionization resonances near the K 
edge, as observed in Figure HI Instead, there is a large number of extremely narrow, yet 



NahaiJ (I2OO5I) 



strong, resonances which need to be fully resolved to account for their area 

uses about 20000 energy points to map out the cross sect ion. For the present calculation 



we use the technique described in 



Witthoeft et al. 



(120091 ) which automatically identifies 



the location of each resonance and adds enough energy points to fully resolve it. Even 
with the effeciency of this approach, over a half million points were required to properly 
resolve all resonance features. To illustrate the strength of the narrow resonances near the 
K edge, we show in the inset of Figure [5] the resonance structure in a narrow energy range. 
The difference between the calculations becomes particularly apparent when compa ring 



convo lved cross sections. Using a Gaussian width of /\E/ E = 10 ^, we find that the 



Nahar 



(I2OO5I ) Ka resonances are about 5 times weaker than in the present results; near the K 
edge, the difference is nearly an order of magnitude. These differences are expected to be 
important for spectral modeling of X-ray plasmas. 

Total photoionization cross sections have been calculated using autostructure 
for Ar-like to Fe-like Ni. The first three ions of this sequence have also been calculated 
using i?-matrix which allows for a comparison of the methods. In Figure [6l we compare 
the present i?-matrix results with two sets of autostructure results using wave 
functions constructed from target or parent orbitals. All data have been convolved with a 
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Gaussian usinp; a width of AE/E = 10 ^. Also shown are the direct photoionization cross 



sections of 



Verner fc YakovlevI (Il995l ). The autostructu re results have 



been s. 



energy (^ 25 eV) so that the K edge is in agreement with 



Palmeri et al. 



lifted in 



fl2008bf ). The 



Verner fc YakovlevI (119951 ) results are unshifted and a disagreement of ~ 50 eV in the 



edge position is observed. The background cross sections above and below the K edge 
from the -R-m atrix and both AUTOSTRU CTURE calculations are in good agreement. The 



background of I Verner fc YakovlevI (Il995l ). on the other hand, is about 10% lower than the 
i?-matrix results just above the K edge, but comes into good agreement by 10 keV. There 
is some disagreement in the position and height of the first resonance (K7) which is due to 
the relaxation effects on the orbitals in the autostructure calculations which are not 
included in i?-matrix. This disagreement becomes smaller for the rest of the resonance 
series. If we were to use orthogonal orbitals in the autostructure calculations, the 
agreement of the resonance positions with i?-matrix would improve, but the discrepancy 
of the resonance heights remains. Since the autostructure calculation using the target 
orbitals gives better agreement with i?-matrix than that using parent orbitals (which is also 
true for Ar-like and K-like Ni), we continue to use target orbitals for the Sc-like to Fe-like 
AUTOSTRUCTURE calculations. Clearly these systems would benefit from future i?-matrix 
calculations, but in the meantime, the autostructure data provide an improvement over 
previous calculations where resonances are neglected entirely. 

By default, autostructure calculates cross sections using the length gauge. This 
begins to fail numerically at high energies above the K edge and the velocity gauge is more 
accurate. We take cross sections calculated using the length-gauge at energies below the 
K edge and the velocity gauge results for energies above the edge. This method yields a 
background cross section in good agreement with the i?-matrix results for Ar-like to Ca-like 
Ni, and is applied to the other systems. 
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Data from this work are to be used in the xstar code (JKallman fc Bautistal 12001 



Bautista &: Kallman 



2OOII ) for modehng photoionized plasmas. To prepare for this 
apphcation, we need to reduce the number of data points from thousands to hundreds. In 
addition to convolution (as described earlier), we remove unneccessary data points from 
each data set in the following manner. Points which lie on a straight line between the two 
adjacent points within 1% are removed. This test is repeated over the entire cross section 
until no more points are removed. Using this process, we can decrease the number of points 
in the cross section data from several thousand to a few hundred, yet still maintain good 
accuracy with linear interpolation. The convolved cross sections are not included with this 
paper, but are available by a request to MCW, xstafo, or the Universal Atomic Database. 
The raw cross sections are available as electronic tables attached to this work. 

The energy range of all datasets in the present work begin at the valence-shell 
threshold. However, as we are only focused on the resonance features near the K edge, we 
have not attempted to fully resolve the resonances in the region of the L edge. 

4. Summary and Conclusions 

Total photoabsorption and photoionization cross sections have been computed for the 
Li-like to Fe-like ions of Ni. Level-resolved, Breit-Pauli i?-matrix calculations were carried 
out for the Li-like to Na-like ions. Term-resolved i?-matrix calculations were performed for 
Mg-like to Ca-like where the mass- velocity and Darwin relativistic corrections are included. 
Distorted wave calculations for Sc-like to Fe-like Ni were performed using autostructure 
where resonances are included using the isolated resonance approximation. 



^http://heasarc. gsfc.nasa.gov /xstar /xstar. html 
^http://heasarc. gsfc.nasa.gov/uadb 
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For all of the i?-niatrix calculations, total photoabsorption and total/partial 
photoionization data are available. Radiation and Auger dampings are included to account 
for radiation losses and the spectator Auger process. No damping is included for the 
AUTOSTRUCTURE calculations so we present only total photoabsorption cross sections which 
are convolved with a Gaussian profile using a width representative of currently deployed 
X-ray detectors. 

We find overa ll good agreement with the background cross sections given by 



Verner fc YakovlevI ( 119951 ). but our positions of the K edge a re typically 40-80 eV higher. 



Naharl ( 120051 ) as to resonance positions. 



For Li-like Ni, we are in good agreement with 

However, we find much stronger resonant enhancement due to the finer energy mesh used 

in the present work. 

All data are provided as on-line tables a ccompanying this article a nd can also be 



obtained through the xstar atomic database ( Bautista &: Kallma: 



20011 ) and the Universal 



Atomic Database. The data sets provided 



and Auger rates reported in 



l ere tog ether with the energy levels and radiative 



Palmeri et al. 



( l2008bl ) will help modelers to carry out detailed 



studies of K spectra of Ni ions. 
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Fig. 1. — Photoabsorption (upper, r^edJ-ancUp^hotoioniza^^ (lower, blue) cross sections of 

Ne-like Ni near the K edge. Color available in the on-line version of the text. 
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Fig. 2. — K/3 reson ance for Ne-lik e^Ms, Si,^, Ar, Qgi, gind Ni. Data from the first five ions 



conies from iWitthoeft et al. 



f>fe^n^'^rie^5V^fe^^ 



(120091 ) ■ The second resonance appearing from Ne-hke S to Ni is 
due to photoexcitation to the Is 2s^ 2p^ 3p ^Pi level. 



-20- 



0.020 



0.015 



1 0.010 
in 

CO 
CO 

o 

6 



0.005 



0.000 




14 16 18 20 

3ar sequence for the present calcu- 



2 4 6 8 10 12 

Fie. 3. — Backeround cr oss sections acrp gs the isQ=ji 

I Electron «.^ount 

lations (filled circles) and IVerner fc YakovlevI ( 119951 ) (open circles). The upper two curves 
(blue) show the cross section above the K edge at 10.9 keV; lower two curves (red) give the 
cross section below the K edge at 6.8 keV. Color available in the on-line version of the text. 
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Fig. 4. — Photoabsorption (top panel) P^J)^tpiqg^^«Ai(fi|^>^V)ttom panel) cross sections 

for Be-like Ni. Red, thin curves show the raw cross section and the blue, thick curves give 
the convolution with a Gaussian having a width of AE/ E = 10^^. Color available in the 
on-line version of the text. 
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Fig. 5. — Total photoionization crp^ section for LijlikgiN i from present, i?- matrix calcu- 

lation (thin, red cu r ve), t he i?- matrix calculation of iNahari ([20051) (thick, green curve), and 



Verner fc Yakovlevi ( 1l995l ) (thin, blue curve with circles ). The Kg resonances are not visible 



in this energy range. The background cross sections of 



Verner fc Yakovlevi ( 119951 ) are in ex- 



cellent agreement with the present work. Inset: close-up of cross sections bet ween 10.08 and 



Nahai] (120051 ) 



10.12 keV; range of cross sections is from 10~^ to 10^. In the inset only, the 

cross sections have been shifted by 8.8 eV to align the resonance positions. Color available 

in the on-line version of the text. 
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Fig. 6. — Photoabsorption cross section fo]^Ca-like,Ni. convolved with a Gaussian having 

a width of AE/ E = 10^^. Based on the height of the resonance feature near 8.45 keV, in 
order from strongest to weakest, the curves are: i?-matrix results (red), AUTOSTRUCTURE 
using target wave functions (green), autqstructure using p arent wave functions (blue). 



and the background cross section of 
the on-line version of the text. 



Verner fc YakovlevI (119951 ) (black). Color available in 



